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AN ANALYTICAL AND EXPERIMENTAL INVESTIGATION 


OF THE MOTION OF A ROTATING SPACE STATION 

By William M. Piland 
Langley Research Center 


SUMMARY 


The dynamics of a manned rotating space-station configuration have been 
investigated both analytically and experimentally. The solutions for the angu- 
lar body motions due to external and internal disturbances were obtained by an 
exact computer method and by an approximate method. The approximate method was 
developed through the linearization of the equations of motion by using small- 
angle theory. In addition, an existing scale model of an inflatable space- 
station configuration was used to obtain experimental results for applied dis- 
turbances on the station. The disturbances considered for this analytical and 
experimental analysis included externally applied torques and static and tran- 
sient products of inertia. 


INTRODUCTION 


In the past, the effects of external and internal disturbances on a 
rotating body in space have been determined by numerical integration of the 
equations of motion with the aid of a digital computer. However, the computer 
solutions can be obtained only for specific configurations and disturbances, 
and they yield little general information. Thus a general closed-form solution 
which would allow a better insight into the effects of the various disturbance 
and configuration parameters would be of considerable value. Furthermore, a 
closed-form method would yield an approximate solution to the equations of 
motion without the aid of computers, even though the closed-form equations 
could also be programed to permit a solution in very short computer time. 

The present report presents the development, results, and use of a closed- 
form solution based on small -angle theory. This solution can be used to pre- 
dict the uncontrolled motion of a spinning vehicle under various disturbances. 
The results obtained from this solution are compared with results obtained from 
a scaled dynamic model which represents a typical space-station conf iguration. 
In addition, these results are checked with the exact solution of the equations 
of motion, as solved by numerical integration with the IBM 709^ electronic data 
processing system. 



SYMBOLS 


a acceleration of moving mass, ft/sec^ 

I moment or product of inertia, slug-ft^ 

m mass of moving tody, slugs 

r radius of moving mass from tody center of gravity, ft 

M external torque, ft -lb 

t time, sec 

t a actual time, sec 

t^ anytime after t s , sec 

t* relative time, sec 

V velocity of moving mass, ft/sec 

X,Y,Z orthogonal body-axis system 

x,y, z, position coordinates of moving mass in body -axis system with origin 
at center of gravity, ft 

v angular velocity of moving mass, radians/sec 

0,9,\|r Euler angles, deg or radians 

to angular velocity, radians/ sec 

m z spin rate, radians/ sec 

Subscripts: 

b body 

f final value 

fs fixed space 

i integer 

max maximum 

o initial value 

s value at the time at which disturbance is removed 
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x,y,z component for X-, Y-, or Z-axis 
xy,xz,yz component for XY-, XZ-, or YZ -plane 
1,2, 3,4 component for mass 1, 2, $, or 4 

A dot over a symbol denotes the derivative with respect to time. 

ANALYSIS 


Presentation of Approximate Solution, Restrictions, 
and Procedure for Use 


The procedure for obtaining an analytical solution for the uncontrolled 
body motion of a rotating spacecraft is developed in the appendix. The effects 
of transient and static product-of -inertia disturbances as well as externally 
applied torques on the body can be analyzed by this procedure. However, in 
order to obtain acceptable accuracy with this solution, certain limitations are 
imposed on the use of the linearized equations of motion. As can be seen in 
the appendix, s mall -angle theory is used in the development of the solution. 
Thus, the modified Euler angles* produced by a disturbance must be less than 
15° • Therefore, in order to apply the approximate solution to a particular 
problem, it must be ascertained that only small angles exist in the problem. 

The maximum angular displacements resulting from a constant product-of-inertia 
disturbance can be approximated by the following equations from reference 1: 


0max 


max 


= tan“I 
= tan “I 


2I yz 
I z - Iy 

2Ixz 

Iz ~ I X 


> 


(i) 


In a manner similar to that of reference 1, the maximum angular displacement 
resulting from a constant applied torque can be approximated by 

, n 2 Mx ^ 

Pma x = tan 


(iz ' Iy) 




( 2 ) 


®max ~ ^ an ^ 


2M-v 


(iz “ Ix) [U 2 


The modified Euler angles used for this analysis are the same as those 
used by Kurzhals in reference 1. As shown in figure 1, the body position is 
defined by three consecutive angles of rotation, \(f, 9, and 0. 
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Equations (l) and (2) are applicable for body configurations with inertia 
ratios Ix/Iz in the interval from 0.5 to 1.0 and zero initial transverse body 
angular velocities and attitude angles. 

The solution presented herein is an example solution and it assumes special 
cases of disturbances. However, most of the basic disturbances that would be 
experienced by manned rotating space stations may be analyzed by this solution. 
For example, assumptions made in the development of the approximate solution 
include a constant velocity of the moving masses, a constant applied disturbing 
torque, and nearly constant spin rates. These assumptions are realistic, but 
are limiting in certain instances. A more general solution which considers all 
possible disturbances can be determined by a method similar to that presented 
in the appendix. 

If a body in space is initially spinning about the axis of maximum inertia 
(the Z-axis) and is disturbed from the steady-spin condition by an externally 
applied moment or an internal mass movement, the motion of the body can be 
determined as shown in the following sections. 

Body motion during the application of disturbance .- For the development of 
the approximate solution, a relative time t* is defined as t* = t a - t Q for 
the interval 0 ^ t a ^ t s (where t s is the time at which the disturbance is 
removed). If the disturbance is applied to the body at zero time, the value of 
t Q is zero and t* is therefore equal to the actual time t a for the interval 
during which the disturbance acts. For this period, the body angular velocities 
about the body X- and Y-axes are 


“x = K 


O - -Meos At* + U o - JjJislaJSi + !Ul + Ji- 

IxA 2 / \ ’ lx* 2 / * lx* lx* 


■* = ■ vfc ) cos Xt * + • 1^)^ + ^ + x 


where 


x* Bot kr 


r* 2 


I, “ Iv 


= - 2 —— 


(tor. 


(3) 


(*) 


(5a) 


V 


Iz ~ I y 

= — = 


* 2 - 


A 1 = “[VV + iyzPz^z + V) + I xz,oV u z^] 


(5b) 

(5c) 

(5d) 


k 


®1 = -Ixz^y^z^ 


(5e) 


A 2 = A X M X + i xz ^ z (fo z + \) - Iy Z ,oVV 

• 2 

B2 - "Iyz^xf°z 

and, for initial spin about the inertial reference axis, 

>> 


°bc,o = 0 

“V,o -0 

_ + (^xz,Q ~ ^-yz^^z) 




“x,o 


• . (iyz,o + Ixz,cP Z )o> z 

°V>o - — + — 


(5f) 

(5g) 


(6) 


The constants A, A]_, Bi, A 2 , B 2 , ux^ 0 , and 0 are dependent only on 

the disturbance applied and the initial characteristics of the body. Therefore, 
the body angular velocities can be solved for directly as a function of time 
for the period during which the disturbance acts. 

Similarly, for the interval 0 ^ t* S t s , the modified Euler reference 
angles can be found from the following equations: 

0 = cos <jo 7 t* + ^2. sin ovt* + — — —(cos At* - cos o> z t*) 

y ° z ^z Z ov 2 -A 2 ' ’ 


D tA 


(o z 2 - A 2 

F 


Ln^tl _ Ew _ cQs A 

V X “z ) Z ' 




avt* - sin 


CJD z t*) 


(7) 
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0 Q \ 

9 = 0 O cos cD z t* + — sin (JD z t* + — 7 — — o( cos _ cos <D z't*) 

^z co z 2 - \ ' 


P2* ( sin At* _ s i n ^ 
cn z 2 - 7\ 2 \ ^ ^z 


+ —^2 (i " cos <0 z t *) 


+ - sin o> z t*^ 


( 8 ) 


T|f = 0> z t* 


(9) 


where the constants are defined as 


. B x A^o z 

°1 - + “**.» - v? ' V? 




„ “> 2^,0 B^ z A X 

D i - ~ir~ - - i^5 + 


Ei =. 2i_ + *25| 

I X A 2 I y A 2 


fl - 


. B ^“z 


IyA 2 


Bp A-I<D„ 

C 2 =^0-^,0 -^ + _ 


/cd^ o . . B x co z A 2 > 

d 2 = -l— t- + >«v,o -r-T-TT 


> 


A 


I x *3 


E 2 


_ B 2 _ Aico z 
' I y A 2 I X A 2 


■p = Bfe 
2 ' “ix^ 2 




(10) 
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and, for initial spin about the inertial reference axis, 

0o = 0 ^ 

9 o - 0 I 

K = “x,o + 'Mo = 0 

9 0 = <*y,o ~ “Z0O = 0 ^ 


(ID 


Body motion after cessation of disturbance . - In the development of the 
solution for the body motion after the disturbance ceases, the same definition 
of the relative time t* was used. Since, by definition, t* = t a - t 0 , then 
at the time the disturbance ceases t Q = t s = t a , and t* in the following 
solution is again initially zero. 

The body angular velocities for any time after time t s are 


where 


^ G / \ 

“x = a bc, s cos to* + s sin + -|(l - cos At*J 


* >2' 


aoy = <jOy^ s cos At* + dby^ s sln fr* + 22.^1 _ cos At*^ 


Gi = 




L xz,f 


Go = - 




•lyz,f 


( 12 ) 

( 15 ) 


w 


and where A x , Ay, and A^ have to be recalculated using the final values of 
the body moments of inertia I x ,f> !y,f> Tz,f‘ These inertias can easily 

be computed if the final position of the moving mass is known. The values of 
d^g and s “ay be found from the following equations by using known 
values of cd^ s and s computed at time t s : 


'^’x, s 
“y,s 


(lz,f ~ I y,f) (D y,^z + lyz^z 2 

J x,f 

(lz,f ~ lx,f) m X,^Z + lxz,f 1> z^ 

x y,f 


( 15 ) 
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The modified Euler angles after time t s are 


<f> = <j) s cos to z t* + < — sin o> z t* 4 — — ^cos “ cos ^z^*) 


a>.* - V 


. ^1^ [ sin At* _ 

<uz 2 - X 2 \ ^ 


sin cD z t‘ 
ov 


-) + ^i 1 ■ oos “ zt *) 


(16) 


®S _ - ... 4-* , ^2 14-* 

s wo uu z u T — sin oj z -g + — 2 “2 


0 = 0o cos o) 7 t* + — sin ovt* + 

“>z o> z * - 7i 


cos At* - cos 


<o z t*) 


+ J g* (s in_M? _ Bln o>,t»\ + KgA . cog t *\ 
oo z 2 - A 2 \ ^ “z / <o z 2 ' ' 


(17) 


i = c^z (t* + t s) 

where 0 S and 0 S are already known and 

■> 

0s = s + ^z^s 

9s = °V,s - ^z0s 1 
-> 

The constants in equations (l6), (17)> and (l8) for <f>, 0, and \|r are 


% = <% s + “z^y, s - 


G^j 


(18) 


(19) 


'Un 


A 2 




(20a) 

(20b) 


Ki = 


G ^°z 


A 2 


TT * G"|0X, 

H 2 = “y,s - “Wx,s + 


(20c) 

(20d) 
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(20e) 



*“sr,« 


s 

X 


K 2 



(20f) 


Experimental Model Description 


An experimental model subjected to typical simulated disturbances was 
analyzed by both the approximate solution and an exact computer solution. The 
rotating model is a l/j-scale model of an inflatable, torus-shaped, space- 
station concept. The full-scale station, shown in figure 2, would be 30 feet 


in diameter and would rotate about the axis of maximum inertia at 6§- rpm to 
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produce a 0.2g environment for the two crew members. The products of inertia 
about the transverse axes of the space station with a mass of 400 slugs would 
be 11 200 slug-ft^ about the X-axis and 12 400 slug-ft^ about the Y-axis. The 
product of inertia about the spin axis would be 1 6 300 slug-ft^. 


The 10-foot-diameter model is constructed of a three-ply neoprene with an 
aluminum hub and has products of inertia about the X-, Y-, and Z-axes of 46, 51> 
and 67 slug-ft^, respectively. The model rotates about the Z-axis at a rate of 
20 rpm. A 4-inch air bearing supports the model and allows ±30° of movement 
about the transverse axes and complete freedom about the spin axis. Disturbance 
torques are applied to the model by the firing of one of four nitrogen jets 
located on the outer periphery of the model, and crew motions or cargo shifts 
are simulated by masses which move in tubes as shown in figure 3- The position 
of these tubes can be changed to yield radial or chordwise mass motions at con- 
stant speeds. For chordwise motion, as shown in figure 3> either of the two 
opposite pairs of masses ^mp and m 2 or and mi^ are moved at the same 

time in such a manner that the crew motion or cargo shift is simulated while 
the center -of -gravity position of the model remains fixed. 


The model is rotated to the desired spin rate by a spin-up motor and shaft 
supported by a boom assembly. The assembly is lifted away from the model after 
the rotational speed has been obtained. The juncture of the spin-up shaft and 
the model is a cone -frustum arrangement which allows the boom to separate from 
the model with little or no tipoff error. 


Instrumentation on the model includes an onboard recorder which traces the 
time history of the body position with respect to inertial space. Onboard 
electrical and nitrogen supplies are also included, and a programing mechanism 
is used to initiate various combinations of disturbance inputs. The position 
of the model is determined by a solar-cell arrangement located on top of the 
hub structure. These solar cells sense the angular deviation of the model with 
respect to a fixed simulated sun on the laboratory ceiling. Angles obtained by 
this arrangement are essentially the small angular displacements of the model 
about the X- and Y-axes. These angular displacements are approximately the 
modified Euler angles, 0 and 9, previously discussed. 
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RESULTS AND DISCUSSION 


The approximate, exact, and experimental results are compared in figures 4 
to 8. These comparisons were made for several disturbances acting on the 
example space station. Disturbances considered included externally applied 
torques and static and transient products of inertia. 


Externally Applied Torque 

Figure 4 shows the motion of the model station due to an externally applied 
torque. The disturbance on the model was produced by firing a jet of 1-pound 
thrust at a moment arm of 5 feet from the model center of gravity for 5 seconds. 
The motion of the model after the torque is removed could simulate the response 
of the full-scale station to the docking impact of a 7000~P°und Gemini vehicle. 
Figure 4(a) shows the time histories of the body angular velocities and toy 

about the transverse body axes and the modified Euler reference angles 0 and 
0. It should be noted that in figure 4(a) is shown a comparison between the 
exact computer solution and the approximate solution, whereas in figure 4(b) is 
shown a comparison between the approximate solution and the experimental results 
obtained from the model. The body angular velocities were not obtained from the 
experimental results since no rate -measuring devices were included in the model 
instrumentation. As can be seen in figure 4, the motions of the model station 
determined by the approximate solution, the exact solution, and the experiments 
are in very good agreement for the applied- torque disturbance, which resulted in 
a maximum angular displacement of 4.75°. 


Static Product of Inertia 


In figure 5 are shown the results of a simulated static or "step" product- 
of -inertia disturbance. Two crewmen are assumed to be located at extreme posi- 
tions in the space station to produce the dynamic unbalance. It should be noted 
that no impulse to simulate the two men instantaneously moving from the center 
of gravity to the extreme positions was imparted to the vehicle. Such a move- 
ment would have been difficult to simulate on the experimental model. However, 
the case presented represents a typical, example of the motion of the station due 
to a crew dynamic unbalance. On the experimental model two 0.111-slug masses 
were located at the following coordinates: x = 3-92 feet, y = 0, and 

z =1.42 feet for mass 1 and x = - 3*92 feet, y = 0, and z = -1.42 feet for 
mass 2. Again, good agreement was obtained from both solutions and the experi- 
mental results; the maximum angular displacement for the unbalance was approxi- 
mately 2.5°. 


Transient Chordwise Products of Inertia 

Figures 6 and 7 illustrate the motion of the station due to transient 
chordwise product-of -inertia disturbances both in the direction of rotation and 
in the direction opposite to the rotation of the station. Movement of the 
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masses in the direction of rotation resulted in a slightly larger maximum angu- 
lar displacement than obtained for movement of the masses in the direction 
opposite to the rotation (approximately I. 7 0 as compared with 1.4°). It was 
assumed that on the full -scale station one 6-slug crewmember moved in a chord- 
wise direction at a rate of 0.6 foot per second for 10 seconds. Shown in fig- 
ures 6 and 7 are "the time histories of the body angular velocities and Euler 
angles of the l/3-scale model for the chordwise movement of two 0.111-slug 

masses at 0.6 foot per second for 3^ seconds. The initial coordinates for the 

location of mass 1 were x 0 = 3*92 feet, y 0 = 0, and z Q = 1.42 feet; those 
for mass 2, x 0 = -3*92 feet, y 0 = 0, and z 0 = -1.42 feet. Each mass moved 
in its respective Z-plane with absolute component velocities of 0.l4l foot per 
second in the x-direction and of O .585 foot per second in the y-direction. 

(See fig. 3.) 

The station motions as calculated by the exact and approximate solutions 
were in very good agreement for both types of chordwise mass movements. Although 
there is a small phase shift between the approximate solution and the experi- 
mental results in the Euler angle plots, the frequencies and amplitudes compare 
very favorably. A possible explanation of these phase shifts is that the inflat- 
able model is not a rigid body. A comparison of the flexible-body motions with a 
rigid-body solution would be expected to result in similar discrepancies. 


Transient Radial Product of Inertia 

In figure 8 are compared the results obtained for the motion of the model 
station due to a transient radial product-of -inertia disturbance. It was 
assumed that on the full-scale station one 6-slug crewman moved radially out- 
ward from the station Z-axis, parallel to and at a constant distance from the 
X-axis. On the model, the two 0.111-slug masses moved parallel to the X-axis 
at a rate of 0.6 foot per second from initial positions xq = 1.62 feet, 
y 0 = 0, and z Q = 1.42 feet for m-j_ and x Q = - 1.62 feet, y Q = 0, and 
z Q = -1.42 feet for m 2 . Both masses moved for 3«1 seconds. Again, as in 
the previous cases of mass-motion disturbances, the frequencies and amplitudes 
of the station motion are in good agreement for both analytical solutions and 
the experimental data obtained. However, the maximum angular deviation experi- 
enced by the station due to the radial mass movement (approximately 1.1°) is 
slightly smaller than that due to the chordwise mass movement. 


CONCLUDING REMARKS 


An approximate solution to the equations of motion for a spinning body in 
space has been developed by using small -angle theory to analyze the effects of 
various disturbances on the body. Results for the approximate solution are com- 
pared with experimental data obtained from a space-station model subjected to 
simulated disturbances and with the results obtained from an exact computer 
solution. These comparisons have shown that extremely good correlation is 
obtainable between the approximate solution and the exact solution. In 
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addition, the results of both analytical solutions have been verified by the 
results of tests performed on the dynamic model. 

The closed-form solution developed in this study may thus be used to 
investigate the effects of both station body characteristics and disturbance 
parameters on the resulting body motion. These investigations will yield rea- 
sonable results for small body angles and rates. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 21 , 19^5* 
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APPENDIX 


DEVELOPMENT OF THE LINEARIZED EQUATIONS OF MOTION 


In order to develop the approximate solution, a detailed analysis was per- 
formed on the exact equations of motion. An example disturbance was used to 
investigate the order of magnitude of the various terms in the equations. The 
particular disturbance considered was the scaled mass of one crewman moving in 
a chordwise direction on the model station as shown in figure 5* The 6-slug 
crewman was assumed to produce a typical transient product -of -inertia distur- 
bance by moving at a constant velocity in the chordwise direction for a total 
of 2 seconds. The scaled mass on the model was divided into two smaller masses, 
mj_ and m 2 , so that the analysis would be simplified by keeping the center-of- 
gravity location stationary. Two other masses, mx and mi^, were present only 

for initial balancing purposes. This type of disturbance analysis also allowed 
easy substantiating of analytical data with results obtained from tests on the 
existing experimental space-station model, the characteristics of which were 
used for this analysis. 

The rotational equations of motion for a spinning body in space, as 
derived in references 1 and 2, are 

Mx = IxPx - Txj4>y - Ixz^z + - txjPy ~ Txz^z 

- <D z ^Iy.Dy - IyzPz - Ixy^bc) + fJi y^Iz ( -°z - IxzPx “ TyzPy) (Ala) 


Wt. • • 

My = ~ -^yz^z “ -^xy^bc ~ ^-yzPz ~ ^xy*bc 

“ ^bc^z^z * ^xz^ ~ ^yz a *y) * m z(^xf°x “ ^xy^V ~ ^xz^z) (Alb) 


M z = I z <i> z - I X z^bc ~ ^-yz^V + ~ ^-xz^x ~ ^yz^y 

- (Wy^I^i>x " Txy°>y - Ixz^z) + “bc^y^y “ %z^z - Ixy°hc) (Ale) 

By definition, the space-station configuration under consideration rotates about 
the axis of maximum inertia (Z-axis) at a spin rate di z . The position of a 

moving mass on the body may be defined at any time by (see ref. 1) 
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x 


y 

z 


12 ^ 
2®x, ot + v x,o t + x o + r cos vt 

|a y ^°t 2 + Vy^ 0 t + y 0 + r sin vt \ 
|a z ,o t2 + V z , 0 t + z Q 


(A2) 


It is not unrealistic to assume that the velocity of the motion of the 
crew is constant, since a crewman walking will accelerate and decelerate almost 
instantaneously. In order to correlate the approximate solution developed 
herein with the results from the experimental model in which linear motion of 
the crew is very closely approximated, the position of the moving mass is 
defined simply as 


x = x 0 + V x Q t 


(A3a) 


y = y Q + v y , 0 t 

and, inasmuch as the motion of the crew is always at the same z-coordinate. 


z = z r 


(A5b) 


If the positions of any number of masses m^ at any time do not affect the 
center-of-gravity position of the overall body, the products of inertia for a 
body are defined as 


•xz - y m i x i 2 i 

(A4a) 

: yz = mjyjZj 

(A4b) 

■xy = y~ m i x iyj 

(Akc) 


In addition, for any number of masses m^, the body moments of inertia about 
the X-, Y-, and Z-axes are defined as 


Ik 
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*x = *x,b +Y1 m i( y i 2 + Z i 2 ) 

x y = + YL mi ( x i 2 + z i 2 ) f 

I z = *z,b + ^ m i( x i 2 + yi 2 ) 


(A5) 


Equations (A5) also apply only if there is no change in the center-of -gravity 
position, as in the case of the experimental model used in the present 
investigation. 


In order to investigate the order of magnitude of the terms of equa- 
tions (Al), the values of the body moments of inertia and products of inertia 
were calculated by substituting the characteristics of the model and of the 
transient product -of -inertia disturbance into equations (A3), (A^), and (A5) • 
The values of these moments of inertia and the derivatives of these inertias 
were then substituted into equations (Al) so that the relative magnitudes of 
the terms of each equation could be analyzed. This analysis was performed at 
the time the masses began to move and again at the time the masses stopped and 
the body disturbance ceased. The effect of the time -varying disturbance on the 
inertia terms were analyzed by using this procedure. Equations (Ala) and (Alb) 
were reduced by assuming that any external torque would be applied only about 
the X- or Y-axis and that the values of and <JOy are much smaller than the 

value of the spin rate <x > z . Eliminating all second-order terms then yielded 


Mx = ^x^x + (iz “ I y) a> z a V " 
My = ly^y " ^z “ lx) “ 


(*xz “ 

> 

(iyZ + IxZ^Z^Z 


(A6) 


and from equation (Ale) the spin speed was found to be approximately constant. 
In addition, the change in the values of the body moments of inertia due to the 
mass movement was negligible. 


Differentiating equations (a 6) and assuming I x , 
stant yields 


ly, I z , and cjd z con- 


Mx = Ixf^x + (iz 
^y = “ (-^z 


ly) (D 7^y 
Ix^z^^x 


(ixz 

C*yz 


+ 


lyz^z) 

ixz^z) 



(A7) 
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and combining equations (a 6) and (A7) yields 

“ AyMy + Ixz^z “ Iyz m z (™z + ^y) “ ^xz^y^z^j 

“y + A^tOy = j=-|^y + >^M X + i'yz^z + ixz^z (^z + ^x) “ lyz^x^z^J 


(A8) 


where 


A* = 


Iz ~ lx 




- Iz " I y- 




A^ - AjjAy 


■0)r, 




(A9) 


Substituting equations (A3)> which represents the positions of the moving 
masses on the model, into equations (A^a) and (A4b) yields 


X xz = mz o( x o + VxjO't) 


Ey Z = mz 0 (y 0 + V J} Q t) 


and the derivatives of these equations are 

i xz = mz 0 V x ^ 0 


Iy Z - mz 0 Vy^ Q 


Ixz = 0 


I yz = 0 


J 


(A10) 


(All) 
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The products of inertia can be written as 


^xz = ^*xz,o + I X z* 


^*yz = l-yz,o + ^yz^ 


(A12) 


Also, the disturbing torques for the experimental model of the present investi- 
gation are not time-varying; therefore M x = My = 0, and, with equations (A12), 
equations (a 8) become 

*\ 


“x + 


~ ~ iyz^z ^°z + V) “ ^xz ,o'V I>z ^ “ ^xzfy^z 2 "^ 

rjiy + A^Oy = ^-JIxMx + ^xz (a z( ta z + A x ) - ~£yz,oty a z^ ~ B yz^xPz^J 


> (A15) 


If a relative time t* is chosen to be t* = t a - t Q for the interval 
0 ^ t a ^ t s (where t s is the time at which the mass stops moving or the 
torque is removed), then, when t Q =0, t* = t a . If this definition of t* 
is used, the body angular velocities and fUy during this interval are 

found from the solution of the second-order differential, equations (A13) 
(assuming a constant spin rate) and are 


OX, = 


(“x, o 


Al 

IxA 2 , 


cos At* + 


B 1 \ sin At* , B l t 




IxA 2 / A 


Z X^ 2 t x a2 


/ (Al4) 


' P’° ■ iSH At * + P'° ■ + ^ + — 


I„A 2 ’ I„A 2 


J 


where 




■ * . » _ i 

Ai = “ j^y^y + -^yz^z^z + ^y) + ^xz^o^y^z J 
®i - -ixzty^z 2 

Ap = A X M X + Ixz^z^z + ^x) “ ^yz^o^x^z 2 
Bp = -Iyz^Uz^ 




(A15) 
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If the tody is initially spinning about the Z-axis only, then 
oo x ^ 0 = ay^ 0 = 0# and from equations (a6) 

/. v 


^,0 

_M X 1 

[-^xz,o ” 

lx 

lx 


_My 1 

(iyz,o + ^xz,cPz) a> z 

— — -f* “ 

% 


(Al6) 


In reference 1, it is shown that the body motion may be related to a set 
of modified Euler angles. The equations involving these angles and the body 
rates are 

0 = + ny tan 0 sin 0 + a> z tan 9 cos 0 

9 = ay cos 0 - a> z sin 0 \ (AIT) 

^ = cn z cos 0 sec 9 + a> sin 0 sec 0 

o ✓ 


Now, if 0 and 0 are assumed to be small, then the sine and tangent of 9 
and 0 can be approximated by the angles themselves and the cosine and secant 
of 9 and 0 can be approximated by unity. By using this small-angle approxi- 
mation and assuming that 90 « 1, 
reduce to 


ay « a> z , and 

a^ « <d z , equations (A17) 

= 0) x + (J0 Z 9 

(A18) 

= ay - m z 0 

(A19) 

= 

(A20) 


The time histories of these modified Euler angles must now be obtained. These 
angles, illustrated in figure 1, define the position of the body with respect 
to a set of fixed-space axes at any time and are necessary to determine the 
body motion completely. Differentiating equation (Al8), solving for 9, and 
substituting into equation (A19) yields 

0 + a> z 2 0 = + uya> z (A2l) 

and in a similar manner the following may be obtained from equations (A19) and 
(Al8) : 
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b" + 0) z ^9 = (By - 0^2 


(A22) 


The equation for \|r can he found by direct integration of equation (A20) . 
After substitution for a^, ay, (%, and ay, equations (A2l) and (A22) are 

solvable differential equations. If the same definition of t* as before, 
during the interval 0 ^ t* ^ t s is used, these solutions yield 

<t) 

6 = 6 cos ax,t* + — sin to 7 t* + i (cos At* - cos u^t*) 

r r o z ox, z m 0 -V? \ z > 


a> z 2 - A 2 ' 


+ fen At* . + 5^( X . cos ^ t .) 

(ii z 2 - A 2 \ A “z / (U Z 2V ' 

+ £i-(ai z t* - sin 0 i z t*) 


(A23) 


9 = 9 0 cos a^ z t* + — sin a> z t* + — ^ o( cos At* - cos ao z t*J 

<J) z a z d - A^ ' 

BgX /sin At* . sin + ft . cos t *\ 
a> z 2 - )>2\ x “z / m^ 2 ' ' 


+ -^r(a> z t* - sin u> z t*\ 


(A24) 


where 


+ = 


C 1 = ^,0 + <W,o 


Bq _ Aga> z 
I X A 2 I y A 2 


^, 0 

D 1 " “A 


A^x,o 


Bgs z Ax 
IyA^ S* 


Ei 


B-j_ A2®z 


(A25) 
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jf _ B 2“z 


^2 = ^y,o “ o 

D 2 = -fe*i£ + 


B 2 + A l a) z 


IyA 2 I X X 2 


B l“ z A 2 


* ' 7 “*'° “ I X A3 " %*/ 


Eo = 


b 2 a 1“z 


IyA 2 I X A 2 


TP _ B l a3 Z 
Jp O 

I x X 2 


A new set of equations are now needed for the station motion after time t s 
when the disturbance is removed. At time t s equations (A13) becomes 


(A26) 



where 




B xz,o + ixz'ts) 

- V Uz2 i - 

~J B xz,f 

X x,f 


0, =-^ 2 | 
2 %,r ' 

B yz,o + B yz t s) 

' 'Iy,f W 

If the time 

t* is again defined as t* = t a 

- t 0 , where now t 0 = t s , then 

the time the 

disturbance is removed, t a = t s 

and the solutions to equa- 


tr ' & ■* 

tions (A26) for and <Uy in the interval t s = t a ^ t ro or 0 ^ t* ^ t 


are: 
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= “x, s cos >** + “x, s + ^(l " cos At*) 


Oiy = oy^ s cos At* + ay^ g (l - cos At*) 




(A27) 


In the equations for o>x and ay, the values of uy^g and oy^ s have previ- 
ously been determined, and from equations (a6). 


“x,s = - 


( x z,f - T y, f )°V, l^Z + Iyz^Z* 


%f 


®JT,B = 


( X Z,f ~ + X xz,f a) z £ 


(A28) 




Substitution of equations (A27) and the derivatives of equations (A27) 
into equations (A21) and (A22) yields two additional second-order differential 
equations which, when solved, are the equations for the modified Euler angles 
0 and 0 after time t s . These solutions are 


0 - 0 f 


cos a> z t* 


+ ijL sin out* + 
(x z z 


(JDr 


— — (cos At* 

,2 - A2 V 


COS (U z 



+ J 1 A ( sin At* 
ay 2 - A 2 \ A 



cos ay 



(A29) 


9 = 0o cos ayt* + — sin out* + ~ (cos At* - cos ayt*) 

S Z ^z Z co 2 _ A 2 V ' 


+ . . J 2^ . + 5s-(i - cos ayt*) <A30) 

a> z 2_A2^ A m z ) ay2 ' Z ' 
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where 


Ghx>„ 

H 1 = “*,8 + “ 2^,8 - 




^yys 


Ki - 


2^2 

X 2 


. GrntD™ 

H 2 = <»y, s " “¥•*, s + -T^ 

A 


= 22 . 




a>zuhc,s 


Gt CO- 
K 2 =-J-5. 

\ 2 

Again at time t* = 0, t a = t s and the values of 9 S and 0 S have already- 
been determined. The values of 0 S and 9 S can be found from equations (Al8) 
and (A19) where, now. 


K = ‘“x, s + “z e 
®s = «V,s - f -°Z0 


-\ 


(A51) 


If the same definition of relative time after time t s is used, the third 
modified Euler angle >|r, from integration of equation (A20), becomes 


♦ = <°z (t* + t s ) 


(a 52 ) 
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Figure 2.- Sketch of inflatable space station. 
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(a) Photograph of experimental model. 
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(b) Schematic drawing of model and moving masses. 
Figure 3.- Characteristics of experimental space-station model. 
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(b) Comparison of experimental results and approximate solution. 
Figure?.- Concluded. 
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(a) Comparison of exact solution and approximate solution. 

Figure 6.- Motion of example station due to a chordwise mass movement opposite to station rotation. 
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Approximate solution 




(a) Comparison of exact solution and approximate solution. 

Figure 7.- Motion of example station due to a chordwise mass movement in the direction of station rotation. 
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(b) Comparison of experimental results with approximate solution. 
Figure 7.- Concluded. 
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(b) Comparison of experimental results with approximate solution. 
Figure 8.- Concluded. 
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